Rheology of ring polymer melts: From linear contaminants to ring/linear blends 
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Ring polymers remain a major challenge to our current understanding of polymer dynamics. 
Experimental results are difficult to interpret because of the uncertainty in the purity and dispersity 
of the sample. Using both equilibrium and non-equilibrium molecular dynamics simulations we have 
systematically investigated the structure, dynamics and rheology of perfectly controlled ring/linear 
polymer blends with chains of such length and flexibility that the number of entanglements is up 
to about 14 per chain, which is comparable to experimental systems examined in the literature. 
The smallest concentration at which linear contaminants increase the zero-shear viscosity of a ring 
polymer melt of these chain lengths by 10% is approximately one-fifth of their overlap concentration. 
When the two architectures are present in equal amounts the viscosity of the blend is approximately 
twice as large as that of the pure linear melt. At this concentration the diffusion coefficient of 
the rings is found to decrease dramatically, while the static and dynamic properties of the linear 
polymers are mostly unaffected. Our results are supported by a primitive path analysis. 



PACS numbers: 83.80.Tc, 83.10.Rs, 83.80. Sg 

While much has been learned about linear and 
branched polymers [1, 2 , a comparable understanding 
of ring or cyclic polymers is lacking. Ring polymers, 
as they do not have free ends, represent the simplest 
model system where reptation is completely suppressed. 
Also, mitochondrial and plasmid DNA are usually cyclic, 
and melts of rings are considered highly relevant model 
systems to understand chromatin folding in the cell nu- 
cleus [21 S] . This makes ring polymers perfect test cases 
for both fundamental polymer and bio physics. 

Early experimental studies on pure ring polymer melts 
gave inconsistent results most likely because the sam- 
ples were contaminated with linear chains [5 - 7 . Also the 
existence of self-knots could not be controlled. More re- 
cently experiments have been conducted [8 based on new 
characterization and purification techniques O [10] . For 
melts of nonconcatenated polystyrene rings with molec- 
ular weight (MW) to entanglement MW ratios of 9.2 and 
11.3, where the entanglement MW is 17500 g/mol, Kap- 
nistos et al. [8 reported that the stress relaxation mod- 
ulus, G(t), follows a power-law decay with no sign of a 
rubbery plateau. The authors used scaling arguments to 
show G{t) ~ t~^l^ ^ a result in agreement with the data 
up to the terminal time. Milner and Newhall [11 intro- 
duced the "diffusion of centrality" concept and mapped 
the ring conformations to annealed tree-like structures 
and found a similar prediction of G(t) ~ Kapnis- 
tos et al. [8 also reported that the smallest concentra- 
tion of linear contaminants that affects the rheology of 
the ring melt is almost two decades below the overlap 
concentration of the linear chains. Despite the synthetic 
effort, the characterization and control of the experimen- 
tal systems including polydispersity, knotting, concate- 
nation and linear contaminants is far from perfect. Be- 
cause of this, computer simulations of optimized models, 



which by now easily reach effective experimental molecu- 
lar weights, are perfect to test concepts for precisely de- 
fined systems under well-controlled conditions. Our own 
recent simulations [12 of a melt of nonconcatenated and 
unknotted ring polymers have shown that G(i) ^ 
with a decreasing from 0.5 to 0.45 with increasing chain 
length. 

Here we employ molecular dynamics (MD) simula- 
tions to study the structure, dynamics and rheology of 
ring/linear polymer blends of equal chain length. We 
consider two lengths of N = 200 and 400 monomers per 
chain. For the model used here the entanglement length 
of a melt of linear polymers is A^e = 28 ± 1 [13 which 
corresponds to N/N^ ^7.1 and 14.3 entanglements per 
chain. For this a bond bending potential {ko = 1.5 e) 
along the chains is introduced, leading to a Kuhn length 
of Ik = 2.79 a [13 , a being the unit of length. is de- 
termined by a primitive path analysis [ISIIH], which is 
known to yield values which properly reproduce rheo- 
logical data [15, 16 . Our systems are perfectly monodis- 
perse, unknotted and nonconcatenated, allowing for a 
rather stringent test of currently discussed concepts. Pre- 
vious simulations of such mixtures have only considered 
short chain lengths and did not measure any rheological 
properties [TTl . While different polymer melts can be 
related to each other by the N/N^ ratio, we note that for 
the present comparison to experiment [8] also the ratio 
of the Kuhn length and the packing length Ij^/p are not 
that different, namely 6.5 for our simulation model (TSj 
and 3.8 for a polystyrene melt [ 19] , 

The topological constraint that a ring must remain 
unknotted and nonconcatenated leads to nontrivial be- 
havior even for the static properties of a melt or con- 
centrated solution of rings. Rings are found to be ap- 
proximately Gaussian at short chain lengths, while for 
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FIG. 1. Universal behavior of {RI{N)) for pure ring poly- 
mer melts. The data were obtained using different simulation 
methods and different models. The reference line with slope 
1/5 corresponds to the Gaussian regime while that with a 
slope of —2/15 corresponds to the collapsed regime. Repre- 
sentative conformations from Ref. ^26^ are shown. 



larger lengths the nonconcatenation dominates the con- 
formational statistics. Gates and Deutsch ^ conjec- 
tured that the exponent in the mean-square gyration ra- 
dius, (i?^) ~ A/"^^, should be less than u = 1/2 and 
greater than 1/3 and used a simple free energy argument 
to arrive at a value of 2/5, which was later supported by 
simulation [21] [22] and experiment [23] for systems with 
less than 13 entanglements per chain. However, for larger 
rings a scaling of {R'^) ~ N'^^^ has been shown [24H26]. 
Altogether we expect a smooth crossover from a Gaus- 
sian regime {u = 1/2) via a regime with = 2/5 for rings 
of length of a few to the "crumpled globule" regime 
{u = 1/3) for rings significantly exceeding TVg. The uni- 
versal scaling behavior of (i?^(7V)) for a pure ring poly- 
mer melt is demonstrated in Fig. [l] using results from 
many different simulations. Only short-chain atomistic 
data for polyethylene [27j deviate from the curve [28] . 
From N/Ne ^ 15 the onset of the collapsed regime is 
clearly observed in agreement with the predictions of Vet- 
torel et al. j24]. 

We present new MD simulations using the same 
semifiexible bead-spring model [29] as in our previous 
work [121 126]. The length, time and energy scales are a, 
r and e, respectively. The production runs were carried 
out using LAMMPS [30] with a time step of 0.01 r and 
an overall monomer density of p = 0.85/cr^. The largest 
simulations ran in parallel on 2048 Blue Gene/P cores. 

Systems studied range from ^Unear = Minear/(Minear + 

brings) = to 1, where M is the number of chains of 
a given architecture. For N = 200 the total number 
of chains ranged from 200-260 while for N = 400 the 
systems were composed of 200-400 chains. The initial 
configuration for each blend system with (/)iinear 

< 0.115 



was created by adding linear chains at random locations 
within an equilibrated ring melt configuration. Ghains 
which most closely matched a Gaussian chain were taken 
from an equilibrated pure linear melt. For the cases 
with (/)iinear ^ 0.25 and 0.5 the appropriate number of 
rings were randomly removed while for the case with 
brings = 10 and Miinear = 250, riugs were taken from 
an equilibrated pure ring melt and inserted into a linear 
melt making sure that the nonconcatenation constraint 
was observed. Because these insertions lead to monomers 
being very nearly overlapping, a short simulation was 
carried out for 100 r while limiting the bead displace- 
ment at every step to 0.001 a. During this short run 
the box size was increased linearly so as to give the cor- 
rect density at the final step. This procedure produces 
non-equilibrated starting configurations. Long MD sim- 
ulations of 4 — 8 X 10^ r were performed to equilibrate 
each system where each architecture moved at least twice 
its root-mean-square gyration radius and in some cases 
more than 20 times this value. 

Results for the mean-square gyration radius for the 
rings and linear chains normalized by their respective 
pure melt values are shown in Fig. |2|a). For the rings 
with TV = 200, {R'^) is found to increase with increas- 
ing linear concentration. At ^unear ^ 0.96, {RD = 
45.3 ± 2.2 cr^, which is 1.5 times larger than the value 
of the pure ring melt. For a Gaussian ring {R'^) = 
Nlkl/ 12 = 45.2 (T^, where / is the average bond length. 
For the rings with ^unear ^0.96 the static structure func- 
tion scales as S{q) ^ q~'^ for 27r/(i?^)^/^ < q < 2ii/lk, 
even though the rings cannot sample the whole confor- 
mational space of a Gaussian ring [31 . For the N = 400 
systems a similar swelling behavior is found for the rings. 
The linear chains are found to be Gaussian for all com- 
binations of N and ^unear- At small values of ^unear the 
rings are partially collapsed as discussed above. As the 
fraction of linear chains increases, the size of the rings 
grows because it is entropically favorable for the linear 
chains to thread the rings. At infinite dilution the non- 
concatenation constraint vanishes and the rings are found 
to be multiply-threaded and nearly Gaussian [TTl |3TJ [32] . 

The diffusion coefficients, which are determined 
by the long-time behavior of the mean-square displace- 
ment of the center-of-mass of the chains, are shown in 
Fig. ^h). The diffusivity of the rings for both val- 
ues of N is found to steadily decrease with increasing 
fraction of linear chains until a dramatic decrease is ob- 
served. With the overlap concentration of linear chains 

being c* = (/>i*nearP = ^/(4/3)7r (i?^ )^^^ , this transition 
corresponds to approximately 0.1 p = 1.5c* for N = 200 
and 0.04p = 0.9c* for TV = 400. For N = 400 the dif- 
fusion coefficient of the rings at ^unear = 0.5 is reduced 
by a factor of about 75 compared to the pure ring melt. 
While the linear chains clearly restrict the motion of the 
rings, the motion of the linear chains for both values of 
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FIG. 2. (a) Mean-square gyration radii, (b) diffusion coeffi- 
cients and (c) zero-sfiear viscosity versus ^linear- Tfie over- 
lap concentration of linear chains, c*, is indicated for the 
two values of N. For the rings with N = 200 and 400, 
{R^)q — 30.8 and 52.9 cr^, respectively, while for the linear 
chains {RI)^ = 88.9 and 180.8 cr^. Note that the horizontal 
axis is interrupted. Lines are drawn as a guide for the eye. 



is largely independent of the blend composition, which 
is consistent with early experimental results [33] , 

Linear chains have free ends and undergo reptation in- 
dependently of whether the surrounding chains are rings 
or linear, and accordingly their diffusion is found to be 
independent of (/)iinear- Rings in a pure melt do not rep- 
tat e like linear chains. As linear chains are added to the 
ring melt, the rings become threaded and the nature of 
their motion changes. A threaded ring can only diffuse 
through the release of threads. For a one-thread situation 
Mills et al. [34] have shown that the diffusion coefficient 
of the ring is ~ ^rTng^iinear- high fractions of lin- 
ear chains the rings become multiply-threaded and their 
diffusion is severely hindered. In this regime the motion 
of a ring monomer is coupled to the motion of surround- 
ing linear chains. This implies Rouse dynamics for the 
ring with a monomer relaxation time governed by the 
reptation relaxation of the linear chains, leading to a re- 



laxation time scaling of ^Hng^unear- "^^^^ argument is 
due to Graessley [3F and predicts D ~ ^rTng^iinear- 

The zero-shear viscosity computed as r^o = G{t)dt 
is shown as a function of (/>iinear 

inFig. [2|c). A striking 
result is the clear indication that the smallest concen- 
tration at which linear contaminants alter the viscosity 
of a ring melt considerably (about 10%) for the chain 
lengths considered here is (/>iinear 

^ 1/100 or c75 with 
a strong increase around c*. This threshold concentra- 
tion is roughly consistent with the change in D for the 
rings. We have confirmed our values of r^o by conduct- 
ing non-equilibrium MD simulations [3F where simple 
steady shear is imposed. For these simulations a Nose- 
Hoover thermostat ^01136] with a relaxation time of 10 r 
was used. Note that the thermal velocity is much larger 
than the largest velocity difference imposed by the shear. 
As shown in Fig. |3]for = 400, when 77(7) is extrapo- 
lated to 7 ^ the agreement with 7^0 is very good j37] . 
Similar agreement is found for N = 200. For both values 
of N the viscosity at (punear = 0.5 is larger than the vis- 
cosity at all other concentrations investigated. For the 
simulated blends with 14.3 entanglements per chain we 

find 7?(^linear = 0.5)/7?o (^linear = 1) ^ 1-8, whcrC T] of 

the blend is taken from the non-equilibrium MD simula- 
tions (cf. Fig. [3| which gives a value that is still increas- 
ing slightly. These findings are in good agreement with 
the experimental results of Roovers [38 who showed for 
ring/linear blends of polybutadiene with approximately 
15.3 entanglements per chain that the maximum in 770 

occurs at (/)linear = 0.6 and 7/o((/>linear ^ 0.5)/7/o ((/>linear = 

1) ^ 2.2. The viscosity results in Fig. [2|c) provide a di- 
rect macroscopic indication of the concentration of linear 
contaminants in experimental samples. As pointed out 
by Kapnistos et al. [Q, the data also suggest how the 
viscosity of a linear melt may be tuned by adding ring 
polymers. 

To quantify the extent of threading, a primitive path 
analysis [131 E] was conducted where the end monomers 
of the linear chains were fixed and the rings were allowed 
to relax freely. This procedure causes the linear chains to 
be pulled taut while the rings shrink towards their cen- 
ter with unthreaded rings collapsing to points. The time 
scale for the primitive path procedure is 10^ r which sat- 
isfies the condition of being equal to or faster than Tg = 
3200 r, the Rouse time of a linear chain of N^. Averaging 
over 10-20 configurations incremented by 2 x 10^ r, with 
N = 200 the percentage of unthreaded rings for ^unear ^ 
0.03, 0.12, 0.25, 0.96 is 80.0, 30.3, 11.7, 0.0%, respec- 
tively. For = 400 with ^unear ^ 0.015, 0.03, 0.12, 0.5 
we find 86.0, 59.0, 7.0, 0.0%, respectively. Fig. ga) 
shows a final configuration for ^unear ^ 0.015 where the 
vast majority of rings are found to be unthreaded. The 
sensitivity of a ring melt to linear contaminants is demon- 
strated by the fact that the viscosity of this system is 
already 1.4 times larger than the pure ring melt value. 
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FIG. 3. Viscosity versus shear rate, 7, for N = 400 obtained 
from non-equilibrium MD simulations [36]. Zero-shear vis- 
cosities obtained from the equilibrium simulations are shown 
on the far left. Note that the horizontal scale is interrupted. 
Inset: Ratio of pure linear to pure ring melt viscosity versus 
number of entanglements per chain for the simulated systems 
and the experimental data of Ref . ^ .39^ . 




FIG. 4. Final configurations from a primitive path analysis 
for N = 400. (a) ^imear = 3/203 ^ 0.015 and ah three hnear 
chains (blue) are shown as well as only the rings which did not 
collapse to points. At this low concentration of linear chains 
on average 86% of the rings are found to be unthreaded, (b) 
^linear = 113/226 = 0.5 and a selected ring (red) is shown 
along with the rings (green) and linear chains (blue) which it 
is either threaded by or entangled with. For clarity all other 
chains are not shown. 



Fig. |4jb) shows a final configuration for ^unear = 0.5 
where a selected ring and the polymers it is entangled 
with are shown. Given the large number of entangle- 
ments at this composition, the dramatic decrease in the 
diffusivity of the rings and the increase in the blend vis- 
cosity in comparison to the pure ring melt value are easily 
understood. 

The present work provides a complete scan of composi- 
tions of two different ring polymer /linear polymer melts 
for dynamical quantities such as viscosity and chain dif- 
fusion. One striking result is that the linear contami- 
nants start significantly affecting the ring melt viscosity 
at a concentration well below their overlap concentration. 
This simulation result is in perfect qualitative agreement 



with the experimental observation of Ref. [8] : according 
to both simulation and experiment, there is clearly an 
effect below the overlap concentration. However, quan- 
titatively we detect the onset of a viscosity change (10% 
increase for rings and linear chains with N/Nq ^ 10) at 
^linear ^ 0.01, while Kapuistos et al. [8] reported a 2- 
fold viscosity increase in comparison to the "pure as cur- 
rently possible rings" at a much smaller concentration of 
</>iinear = 0.0007. To providc an intuitive picture of these 
concentrations one can estimate what would be the typ- 
ical distances between chains. For c^)iinear = 0.0007 the 
typical distance between linear chains {p/N)~^^^ would 
be about 66 a for N = 200 and 83 a for TV = 400. The 

1/2 

diameter {2{Rl) ' ) of the rings is about 11 a and 15 a 
and of the linear chains about 19 a and 27 a, respec- 
tively. Thus two linear chains would be separated on 
average by about 4-5 ring diameters for N = 200, or by 
4 ring diameters for N = 400. And, importantly, these 
rings would not be entangled since they are nonconcate- 
nated and have no free ends. At (punear = 0.01, where 
our data indicate a 10% viscosity increase, distances and 
chain extensions are all rather similar. 

While the two works differ with respect to the onset 
concentration, fair agreement is found for the ratio of the 
pure linear melt viscosity to that of the (almost) pure ring 
melt as shown in the inset of Fig. [s] [39] . Additionally, 
the simulation and experimental results for the depen- 
dence of ?70, linear /^O, rings OU N/Nq are COUSistCUt UOt Ouly 

with one another, but also with the theoretical framework 
of Ref. [8] and our previous result [12] which suggest a 
power law dependence with power close to 2. 
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